Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species’ range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals

Evolutionary and demographic correlates of Pleistocene coastline changes in the Sicilian wall lizard Podarcis wagleriana

Aim Emergence of coastal lowlands during Pleistocene ice ages might have provided conditions for glacial expansions (demographic and spatial), rather than contraction, of coastal populations of temperate species. Here, we tested these predictions in the insular endemic Sicilian wall lizard Podarcis wagleriana. Location Sicily and neighbouring islands. Methods We sampled 179 individuals from 45 localities across the whole range of P. wagleriana. We investigated demographic and spatial variations through time using Bayesian coalescent models (Bayesian phylogeographic reconstruction, Extended Bayesian Skyline plots, Isolation‐with‐migration models) based on multilocus DNA sequence data. We used species distribution modelling to reconstruct present and past habitat suitability. Results We found two main lineages distributed in the east and west portions of the current species range and a third lineage restricted to a small area in the north of Sicily. Multiple lines of evidence from palaeogeographic (shorelines), palaeoclimatic (species distribution models), and multilocus genetic data (demographic and spatial Bayesian reconstructions) indicate that these lineages originated in distinct refugia, located in the north‐western and south‐eastern coastal lowlands, during Middle Pleistocene interglacial phases, and came into secondary contact following demographic and spatial expansions during the last glacial phase. Main conclusions This scenario of interglacial contraction and glacial expansion is in sharp contrast with patterns commonly observed in temperate species on the continent but parallels recent findings on other Mediterranean island endemics. Such a reverse expansion–contraction (EC) dynamic has been likely associated with glacial increases of climatically suitable coastal lowlands, suggesting this might be a general pattern in Mediterranean island species and also in other coastal regions strongly affected by glacial marine regressions during glacial episodes. This study provides explicit predictions and some methodological recommendations for testing the reverse EC model in other region and taxa

Climate change promotes hybridisation between deeply divergent species

Rare hybridisations between deeply divergent animal species have been reported for decades in a wide range of taxa, but have often remained unexplained, mainly considered chance events and reported as anecdotal. Here, we combine field observations with long-term data concerning natural hybridisations, climate, land-use, and field-validated species distribution models for two deeply divergent and naturally sympatric toad species in Europe (Bufo bufo and Bufotes viridis species groups). We show that climate warming and seasonal extreme temperatures are conspiring to set the scene for these maladaptive hybridisations, by differentially affecting life-history traits of both species. Our results identify and provide evidence of an ultimate cause for such events, and reveal that the potential influence of climate change on interspecific hybridisations goes far beyond closely related species. Furthermore, climate projections suggest that the chances for these events will steadily increase in the near future

How complex is the Bufo bufo species group?

Species delineation remains one of the most challenging tasks in the study of biodiversity, mostly owing to the application of different species concepts, which results in contrasting taxonomic arrangements. This has important practical consequences, since species are basic units in fields like ecology and conservation biology. We here review molecular genetic evidence relevant to the systematics of toads in the Bufo bufo species group (Anura, Bufonidae). Two studies recently published in this journal (Recuero et al., MPE 62: 71-86 and García-Porta et al., MPE 63: 113-130) addressed this issue but reached opposing conclusions on the taxonomy of the group (four versus two species). In particular, allozyme data in the latter paper were interpreted as evidence for hybridization across species (between B. bufo-B. spinosus and B. bufo-B. verrucosissimus). We tested claims for hybridization through re-analysis of allozyme data for individuals instead of populations, to be able to distinguish between sympatry with and without admixture, and found no evidence of hybridization across taxa. We propose alternative explanations for the observed patterns that García-Porta et al. (2012) failed to consider. In the absence of unequivocal evidence for hybridization and introgression, we reject the proposal to downgrade Bufo spinosus and Bufo verrucosissimus to the subspecies level.IMS is a "Ramón y Cajal" postdoctoral fellow supported by the Spanish Ministerio de Ciencia e Innovación and the Universidad de Castilla la Mancha. Partial funds for this work were provided by grants CGL2008-04271-C02-01/BOS and CGL2011-28300 (Ministerio de Ciencia e Innovación, Ministerio de Economía y Competitividad, FEDER) and PPII10-0097-4200 (Junta de Comunidades de Castilla la Mancha) to IMS.Peer Reviewe

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae)

et al.New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490. base. pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.IMS is a “Ramón y Cajal” postdoctoral fellow supported by the Spanish Ministerio de Ciencia e Innovación and the Universidad de Castilla la Mancha. The Synthesis Project (http://www.synthesys.info/) of the European Union partially supported this study (Ref.: HU-TAF-181). Partial funds were also provided by the Spanish Ministerio de Ciencia e Innovación (Ref.: CGL2008-04271-C02-01/ BOS) and Junta de Comunidades de Castilla la Mancha (Ref.: PPII10-0097-4200) to IMS, Hungarian Scientific Research Fund (OTKA K77841) to JV and Ministry of Education and Science of Republic of Serbia (Grant No. 173025) to JCI.Peer Reviewe

Drilling down hotspots of infraspecific diversity to bring them into on-ground conservation of threatened species

Unprecedented rates of biodiversity loss raise the urgency for preserving species ability to cope with ongoing global changes. An approach in this direction is to target intra-specific hotspots of genetic diversity as conservation priorities. However, these hotspots are often identified by sampling at a spatial resolution too coarse to be useful in practical management of threatened species, hindering the long-appealed dialog between conservation stakeholders and conservation genetic researchers. Here, we investigated the spatial and temporal variation in species presence, genetic diversity, as well as potential risk factors, within a previously identified hotspot of genetic diversity for the endangered Apennine yellow-bellied toad Bombina pachypus. Our results show that this hotspot is neither a geographically homogeneous nor a temporally stable unit. Over a time-window spanning 10–40 years since previous assessments, B. pachypus populations declined in large portions of their hotspot, and their genetic diversity levels decreased. Considering the demographic trend, genetic and epidemiological data, and models of current and future climatic suitability, populations at the extreme south of the hotspot area still qualify for urgent in-situ conservation actions, whereas northern populations would be better managed through a mix of in-situ and ex-situ actions. Our results emphasize that identifying hotspots of genetic diversity, albeit an essential step, does not suffice to warrant on-ground conservation of threatened species. Hotspots should be analyzed at finer geographic and temporal scales, to provide conservation stakeholders with key knowledge to best define conservation priorities, and to optimize resource allocation to alternative management practices

Genomic Evidence for Cryptic Speciation in Tree Frogs From the Apennine Peninsula, With Description of Hyla perrini sp. nov

Despite increasing appreciation of the speciation continuum, delimiting and describing new species is a major yet necessary challenge of modern phylogeography to help optimize conservation efforts. In amphibians, the lack of phenotypic differences between closely-related taxa, their complex, sometimes unresolved phylogenetic relationships, and their potential to hybridize all act to blur taxonomic boundaries. Here we implement a multi-disciplinary approach to evaluate the nature of two deeply-diverged mitochondrial lineages previously documented in Italian tree frogs (Hyla intermedia s. l.), distributed north and south of the Northern Apennine Mountains. Based on evidence from mitochondrial phylogenetics, nuclear phylogenomics, hybrid zone population genomics, niche modeling analyses, and biometric assessments, we propose that these lineages be considered distinct, cryptic species. Both mitochondrial and nuclear data affirm that they belong to two monophyletic clades of Pliocene divergence (~3.5 My), only admixing over a relatively narrow contact zone restricted to the southeast of the Po Plain (50–100 km). These characteristics are comparable to similarly-studied parapatric amphibians bearing a specific status. Inferred from their current geographic distribution, the two Italian tree frogs feature distinct ecological niches (<15% of niche overlap), raising questions regarding potential adaptive components contributing to their incipient speciation. However, we found no diagnostic morphological and bioacoustic differences between them. This system illustrates the speciation continuum of Western-Palearctic tree frogs and identifies additional cryptic lineages of similar divergence to be treated as separate species (H. cf. meridionalis). We recommend combined approaches using genomic data as applied here for the future taxonomic assessment of cryptic diversity in alloparapatric radiations of terrestrial vertebrates, especially in controversial taxa. Finally, we formally described the northern Italian tree frogs as a new species, Hyla perrini sp. nov

An extended mtDNA phylogeography for the alpine newt illuminates the provenance of introduced populations

Many herpetofauna species have been introduced outside of their native range. MtDNA barcoding is regularly used to determine the provenance of such populations. The alpine newt has been introduced across the Netherlands, the United Kingdom and Ireland. However, geographical mtDNA structure across the natural range of the alpine newt is still incompletely understood and certain regions are severely undersampled. We collect mtDNA sequence data of over seven hundred individuals, from both the native and the introduced range. The main new insights from our extended mtDNA phylogeography are that 1) haplotypes from Spain do not form a reciprocally monophyletic clade, but are nested inside the mtDNA clade that covers western and eastern Europe; and 2) haplotypes from the northwest Balkans form a monophyletic clade together with those from the Southern Carpathians and Apuseni Mountains. We also home in on the regions where the distinct mtDNA clades meet in nature. We show that four out of the seven distinct mtDNA clades that comprise the alpine newt are implicated in the introductions in the Netherlands, United Kingdom and Ireland. In several introduced localities, two distinct mtDNA clades co-occur. As these mtDNA clades presumably represent cryptic species, we urge that the extent of genetic admixture between them is assessed from genome-wide nuclear DNA markers. We mobilized a large number of citizen scientists in this project to support the collection of DNA samples by skin swabbing and underscore the effectiveness of this sampling technique for mtDNA barcoding

Multilocus Bayesian Estimates of Intra-Oceanic Genetic Differentiation, Connectivity, and Admixture in Atlantic Swordfish (Xiphias gladius L.)

Versión del editor